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ABSTRACT 


A  two-step  combustion  process,  consisting  of  monopropellant  (i.e.  one- 
reactant)  burning  followed  by  bipropellant  (i.e.  two-reactant)  burning,  is 
considered.  The  original  reactant,  namely  vapor  from  a  liquid  droplet,  decom¬ 
poses  into  a  fuel  that  combines  with  oxidant  in  the  surrounding  atmosphere . 
Complete  responses  M(D^,D2)  of  the  evaporation  rate  M  of  the  droplet  to 
the  DamkOhler  numbers  ,D2  of  the  two  reactions  are  determined  in  the  limit 
of  large  activation  energies.  Conditions  under  which  the  response  is  ntono- 
tonic  or  multi-valued  (thereby  exhibiting  auto-ignition  and  auto-extinction) 
are  identified.  Previous  conjectures,  based  on  DamkOhler-number  asymptotics, 
are  found  to  be  mostly  correct. _ 
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SIGNIFICANCE  AND  EXPLANATION 


Asymptotic  Methods  are  responsible  for  remarkable  advances  in  combustion 
theory  over  the  last  decade,  however  they  have  been  largely  limited  to  single- 
step  reactions.  This  paper  uses  such  methods  to  investigate  an  example  of 
multiple-step  reactions,  which  are  of  greater  practical  interest. 

The  following  two-step  process  is  considered.  Vapor  given  off  by 
evaporation  of  a  liquid  droplet  decomposes  into  a  fuel.  (This  decomposition 
constitutes  the  first  step  of  the  combustion. )  Bie  fuel  produced  then  reacts 
with  the  surrounding  oxidizing  atmosphere  in  the  second  step.  The  burning  of 
a  hydrazine  droplet  in  oxygen  is  an  exa^le  of  such  a  reaction. 

The  goal  of  this  investigation  is  to  determine  the  evaporation  rate  since 
it  measures  the  rate  of  consumption  of  the  reactants,  a  quantity  of  primary 
interest  in  the  applications  of  combustion  theory.  This  is  accomplished  by 
using  the  method  of  matched  asymptotic  expansions  based  on  the  physically 


realistic  limit  of  large  activation  energies.  Evaporation  rates  which  are 


proportional  to  the  droplet  surface  area  or  to  its  radius  are  found.  In 
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addition,  conditions  for  ignition  and  extinction  are  identified.  DTIC  TAB 
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DROPLET  DECOMPOSITION  IN  A  REACTIVE  ATMOSPHERE i 
COMPLETE  RESPONSES  FOR  LARGE  ACTIVATION  ENERGIES 

•  •• 

H.  V.  McConnaughey  and  G.  S.  s.  Ludford 

1.  Introduction 

The  central  question  of  coatbustion  theory  is  to  determine  how  fast  the  reactants  are 
con aimed.  When  there  is  a  single-step  reaction  this  amounts  to  determining  the  burning 
rate  M  as  a  function  of  the  so-called  DamkShler  number  Di  other  parameters  may  be 
involved ,  but  they  are  supposed  to  be  held  fixed.  For  example,  consider  the  bipropellant 
(i.e.  two-reactant)  burning  of  a  fuel  droplet,  in  which  the  liquid  in  a  spherical  droplet 
evaporates  at  its  surface  and  tha  vapor  is  oxidised  by  the  ambient  atmosphere  in  a 
concentric  reaction  tone.  The  DamkBhler  number  depends  notably  on  the  rate  of  chemical 
reaction,  the  ambient  pressure,  and  the  radius  of  the  droplet!  any  of  these  can  be  used  to 
vary  D.  The  response  M(D)  is  to  be  found  for  various  values  of  three  additional 
parameters!  the  droplet  temperature  Ta,  the  ambient  temperature  T^,  and  the  latent 
heat  of  evaporation  L.  Tha  shape  of  the  response  curve  depends  on  the  (fixed)  values  of 

T  -  T  and  L. 

•  a 

In  general,  the  analytical  determination  of  M(D)  is  prevented  by  the  nonlinearity  of 
the  governing  equations,  which  persists  in  the  (Arrhenius)  reaction  term  even  when  simple 
geometries  (such  as  that  of  the  droplet  problem)  are  adopted.  The  Arrhenius  nonlinearity 
has  bean  overcame  by  asymptotic  methods,  in  which  the  activation  energy  0  tends  to 
infinity;  this  is  a  physically  realistic  limit  since  many  reactions  of  interest  in 
combustion  do  have  large  activation  energies.  For  the  fuel  droplet  this  leads  to 
mono  tonic,  S- shaped,  and  (for  the  practically  unimportant  case  T—  <  Tg )  C-shaped  response 
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curves,  depending  on  the  values  of  T—  -  and  L.  Similar  results  are  found  for  the 
no nopropellant  (i.a.  a ingle-reactant)  burning  of  other  droplets,  in  which  the  vapor 
decomposes  into  products  in  a  concentric  reaction  zone.  (See  Buckmaster  and  Ludford  [1, 
pp.  100  and  127].) 

Asyag>totic  methods  are  responsible  for  raawrkable  advances  in  coa&ustion  theory  during 
the  last  decade,  even  though  they  were  largely  limited  to  single-step  reactions.  A  recent 
account  of  some  of  these  successes  has  been  given  by  Buckmaster  and  Ludford  [2].  Neverthe¬ 
less,  the  need  for  a  sxjra  thorough  treatment  of  multiple-step  reactions  has  been  recognized 
for  some  time,  and  isolated  efforta  are  now  being  replaced  by  broad  action. 

The  simplest  problem  of  cca&ustion  theory  is  the  steady,  unbounded  plane  fleam,  and  it 
is  natural  to  begin  an  investigation  of  multiple-step  reactions  there.  That  restricts  the 
discussion  to  premised  flames,  i.e.  the  combustion  of  reactants  that  are  already  mixed, 
because  plane  diffusion  flames,  in  which  the  reactants  are  originally  separate  and  must  mix 
by  diffusion,  only  exist  when  limited  in  extent.  In  order  to  extend  the  existing 
discussion  to  diffusion  flames,  we  have  therefore  adopted  a  spherical  geometry  which, 
because  it  is  still  one-dimensional,  introduces  a  minimus  of  complications. 

The  multi-step  reaction  that  we  have  chosen  to  examine  is  a  combination  of  the  mono- 
propellant  and  bipropellant  reactions  mentioned  above  that  also  has  practical  interest. 
Vapor  from  a  droplet  decomposes  at  a  premixed  flame  into  a  fuel  (and  possibly  other 
products)  that  is  oxidised  by  the  ambient  atmosphere  at  a  diffusion  flame  (Figure  1).  The 
complete  response  of  the  burning  rate  to  the  two  Da mk Shier  numbers  D^D.,  io  determined  as 
the  two  activation  energies  91 ,8^  tend  to  infinity.  (An  atmosphere  hotter  than  the 
droplet,  i.e.  ?m  >  T(,  is  assumed  since  the  reverse  is  not  of  great  practical  interest.) 

A  popular  experimental  study  of  this  type  of  hybrid  combustion  investigates  the 
burning  of  a  hydrazine  droplet  in  an  oxidizing  environment.  Variation  of  the  evaporation 
rate  with  droplet  diameter  or  with  pressure  are  typically  measured  (Dykema  and  Greene  [3], 
Lawver  [4],  Dosser  and  Peskin  (5],  Allison  and  Faath  [6]).  Hybrid  droplet  combustion  is 
also  mentioned  by  Williams  [7,  pp.  246-7],  who  points  out  the  need  for  theoretical 
investigation  of  the  process.  Our  results  can  hardly  be  compared  quantitatively  with  such 


experiments,  but  their  qualitative  feature*  suggest  that  acre  extensive  experiments  should 


be  made . 

Theoretical  work  on  the  problem  has  been  done  by  Fendell  [8]  and  by  Buckmaster,  Kapila 
and  Ludford  [9].  Both  these  investigations  are  based  on  Damkfihler-number  asymptotics,  i.e. 
large  or  seall  and  large  or  small  D 2»  hence  they  can  described  at  most  the  "corners" 

of  the  complete  response  surface  M(D.j  ,D2 ) .  Fendell  identifies  the  resulting  burning-rate 
formulas,  while  Bucknaster  et  al.  clarify  and  add  to  Fendell* s  results,  suggesting  how  the 
valid  corner  values  are  connected  by  the  "edges"  of  the  response  surface  M(D^,d2).  The 
latter  authors  acknowledge  the  tentativeness  of  their  findings  and  the  need  for  analysis 
with  finite  D.|  and  D2,  recognising  that  "such  analyses  can  be  carried  out  when  the 
activation  energies  6^  and  8^  are  large,  but  such  a  description  is  quite 
ccs^licated. . 

Details  of  the  complication  have  been  excluded  from  our  account  as  far  as  possible. 

Of  the  18  types  of  asymptotic  solutions  we  describe  just  3  typical  ones,  relegating  the 
supporting  analysis  (in  the  main)  to  impressionistic  appendices.  However,  the  tables  and 
figures  do  give  a  complete  picture  of  the  results  and  further  information  can  be  found  in 
McConnaughey's  thesis  [10].  Finally,  the  results  are  shown  to  agree  with  earlier,  well 
established  findings,  and  are  compared  and  contrasted  with  the  conjectures  of  Buckmaster  et 
al.  [9]. 


2.  Formulation 

The  problem  formulation  and  notation  of  Buckmaster  et  al.  [9]  is  adopted.  The 
governing  equations  in  dimensionless  form  are 

L(Tm>  “  D1*V  t(V  “  _Dl“l  +  D2W2'  l(Yo>  "  W  MT)  "  "‘Vi  ”  QD2“2  '  (2,1) 

where 


1  =  d2/dr2  ♦  [(2r-M)/r2]d/dr,  »1  -  Y^  expt-Bj/T),  ci>2  -  YfYQ  exp(-82/T>  , 

and  the  other  quantities  are  defined  in  the  list  of  symbols .  The  boundary  conditions  are 


1: 


dY 

a 

dr 


A\ 


dY 


dT 


MY  ,  7s  “  MY  ,  j 
f  dr  o  dr 


ML,  T 


(2.2) 


(2.  3) 


«■  "  *,  -  rf  -  o,  ro  -  V'  T  "  T-  "  T.  • 

A  derivation  of  this  modal  ia  sketched  by  McConnaughay  [10],  based  on  a  more  general 
development  given  by  Bucknaster  and  Ladford  [1]. 

equations  (2.1)  can  be  combined  and  integrated  subject  to  (2.2)  and  (2.3)  to  yield  the 


Shvab-Xeldovich  relations 


♦  Tf  -  Yo  -  1  -  d  ♦  Y^H 


Y  +  QT  +T-I  -  L  +  1  ♦  (T  -  T  ♦  L  -  1  ♦  QY  >e”^  ,  (2.5) 

■  O  I  •  I  o* 

which  are  vaiid  everywhere.  Thus  only  two  of  equations  (2.1)  are  needed.  We  select 
equations  (2.1a)  and  (2. Id)  and  the  corresponding  boundary  conditions  to  form  a  fourth- 
order  system  subject  to  five  boundary  conditions,  the  unknown  constant  M  is  thereby 
determined  as  a  function  of  the  various  parameters  which  appear. 

The  objective  of  this  work  is  to  describe  the  behavior  of  M(D^,D2)  for  all  values  of 
its  arguments,  by  considering  the  independent  limits  8^  ♦  “  and  0^  ♦  •». 


In  accordance  with  the  usual  strategy  of  activation-energy  asymptotics  (Buckmaster  and 
Ladford  [1]f  Kepi la  and  Ladford  [11]),  the  DamkQhler  numbers  are  written 

expdjA^),  i  «  1  or  2  ,  (3.1) 

where  is  a  positive  paraamter  and  is  at  most  algebraic  in  8^.  The  ith  reaction 

term  then  contains  the  expression  D*  axp( 8^ A^-H^/'T)  and  the  limit  8^  ♦  •  confines  the 
associated  chemical  activity  to  a  thin  layer  either  at  temperature  Ti#  «  Ti  with 
thickness  of  order  T^/8^  or  at  temperature  T^#  >  with  thickness  of  order 
exp[c(8i/Ti<,-8i/Ti)] ,  where  c  “  1/2  or  1/3.  The  reactions  are  negligible  outside  of 
these  flame  sheets,  hence  the  majority  of  the  cos&ustion  field  is  described  by  the  linear 


system 


L(Y  )  -  L(T)  -  0  . 


Matched  asymptotic  expansions  are  then  used  to  describe  the  solution  inside  and 
outside  the  flasw  cones. 
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3.1 


The  solution  obtained  by  this  method  depends  on  the  assumed  flame  oonf  iguratioa  and 
flame  types.  With  the  location  of  the  decoatpoaition  reaction  denoted  by  r  1  and  with  the 
bipropellant  flame  at  r2,  It  ie  found  that  all  possible  configurations  for  which 
1  <  r1  <  r^  <  •  admit  a  solution.  Whether  or  not  Yf(ra>  and  ro(r2)  vanish  as 

6^  ♦  determines  the  nature  of  the  flames,  which  ia  often  not  unique  for  a  given  flame 
arrangeawnt .  Bence,  for  each  configuration,  it  is  necessary  to  consider  separately  the 
various  types  of  flames  possible,  i.e.,  partial  vs.  complete  decomposition  at  r,  and  a 
complete-burning  (or  so-called  Burke-SchUMnn)  diffusion  flame  va.  a  partial  burning  flame 
at  r2.  When  all  possible  ccafeinationa  of  flame  locations  and  types  are  investigated,  the 
cases  shown  in  Table  1  are  found  to  admit  solutions  of  interest.  (Other  conceivable  cease, 
see  McConnaughey  (10],  are  omitted)  they  correspond  to  special  cases  whose  treatment  has 
been  omitted  for  simplicity.)  The  first  thirteen  categories  represent  separated  flames  for 
which  M  -  0 ( 1 )  i  the  next  four  correspond  to  merged  flames  with  M  -  0(1);  the  last  is 
the  only  configuration  for  N  >>  1.  The  asymptotic  analysis  of  these  different 
possibilities  is  illustrated  in  the  following  sections  and  associated  appendices. 

3.2.  Separated  names  with  H  ■  0(1) 

Solution  category  6,  which  is  illustrated  in  Figure  2,  is  considers.*  here.  The 
reactionless  system  (3.2)  is  valid  for  1  <  r  <  r)(  r,  <  r  <  r2  and  r2  <  r  <  •»  the 
leading  terns  of  the  corresponding  asymptotic  expansions  are 

1  <  r  <  r1«  T ^  «  1  -  e  1  ,  T  ■  T  -  L  ♦  LeMe~M/r  ,  (3.3) 

r,  <  r  <  r,»  T  -  0,  T  -  T  -  L  ♦  1  ♦  (T  -  T  ♦  L  -  1  ♦  QY  >e_H/r  ,  (3.4) 

1  4  m  1  ■  • 

r,  <  r  <  -  J  Y»0,  T-T-L  +  1+  Q  +  (T  -  T  ♦  L  -  1  -  Q)a~ti/x  , 

4  ■  8  •  fl 

where  conditions  and  relations  (2.2)  -  (2.5)  have  been  invoked  as  have  Ya(r1)  »  o(1), 

Yf ( r2 )  "  o(1),  Y0(r2)  -  o(1)  and  continuity  to  leading  order  of  VB,  Yf  and  Yq  at  r, 
and  r2> 
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tkble  1 


solution  of  interest 


Solution 

Category 

Description 

1 

■  1  «  r2  «  »i  T^(rg)  ■  o(l)  i  YQ(r2)  *  o(l) 

2 

»  1  <  ra  «  »|  T^(rg)  «  o(l)i  70(r2)  ■  0(1) 

3 

»  1  «  r2  <  •»  t^(r2)  ■  0(1) »  To(r2)  »  o(l) 

k 

Tx  •  1  <  r2  <  «j  If (rg)  ■  0(1) %  7o(r2)  »  0(1) 

5 

rl“1,r2-" 

6 

1  *  rl  K  *2  *  “*  T«^rl^  *  Xt^T2^  “  °^*  Y</r2* 

»  o(l) 

T 

1  <  rx  <  r2<  •!  7f(rg)  »  o(l)s  Y0(r2)  -  0(1) 

8 

1  <  rx  <  r2  <  -i  TMr^  .»  o(l)»  If(r2)  »  0(1)  i  yo<r2) 

-o(l) 

9 

1  <  rl  <  r2  *  "*  Tf*r2^  *  7^(r2)  *  0(1) 

10 

1  *  rl  <  t2  *  *!  T»^rl^  *  ° M 

11 

1  <  rx  <  rg  <  •»  Y^Tj)  »  0(1)  i  Tf(r2)  «  o(lh  YQ(r2) 

-  o(l) 

12 

1  <  rl  *  r2  *  "l  Tm^rl^  "  *f(r2)  *  °^;  To^r2^ 

■  0(1) 

13 

1  <  rx  <  rg  ■  Y^)  ■  0(1) 

Ik 

1  *  rl  *  r2 

15 

1  <  r^  ■  r2  <  •»  Y^r^)  ■  o(l)»  Yf(rg)  ■  o(l)v  YQ(r2) 

«  0(1) 

18 

1  <  rx  »  r2  <  •»  T^r^)  ■  0(1)  i  Yf(r2)  ■  o(lh  Y^(r2) 

•  0(1) 

17 

rA  -  r2  -  - 

18 

N  »  li  r,  ■  1»  r.  •  • 

r 


Sine*  Ya(r)  ¥  0  for  r  <  ty,  the  temperature  must  be  lese  than  T1  in  that  region, 
for  otherwise  the  reaction  term  exp(8j/Yj-6^/T)  would  not  be  negligible  there.  In 

order  for  the  reaction  to  become  important  at  r^  then,  we  auet  have  T(r^)  -  T1  >  Tg. 
Given  this  and  continuity  to  leading  order  of  T  at  r,,  and  r2,  the  solution  shows  that 


H/r,  T,yra*L-HQYe 


T,-T  +L-1 

1  S 


,T1“T«+L  M/r,  */r2 

i — r— ]•  ,  • 


-  i  +  r 


The  outer  expansions  alone  have  therefore  determined  M,  r1  and  r2  as  functions  of 
various  parameters.  In  order  that  these  results  be  consistent  with  the  hypothesis 
1  <  r1  <  r2,  the  following  conditions  oust  hold: 

t  >  0»  T^-L+1  <  Tt  <  T(r2)  or  T(r2>  <  T<  <  T#-L*1  ,  (3.5) 

where  T(r,)  -  [T  +Y  (T  -L+1+Q))/(1+Y  >  S  T  . 

It  is  also  necessary  to  examine  the  structure  of  the  two  flames  in  order  to  discover 
what  additional  conditions  restrict  the  validity  of  the  above  results.  Such  restrictions 
amount  to  existence  conditions  for  solution  of  the  inner  problems  at  r.  and  r2,  where 
the  expansions  must  match  the  outer  ones.  The  associated  analysis  is  sketched  in  Appendix 
1  and  adds  the  requirements 

T2  <  TB8;  T1  >  VL+  2  "h,n  Tl  >  TBS  * 

This  coagtletes  the  investigation  of  category  6.  In  a  manner  similar  to  (.hat  just 

demonstrated,  the  solution  in  each  of  the  separated- f 1 asm  categories  1-13  along  with 

restrictions  on  its  validity  are  determined.  The  burning  rates  are  listed  in  Table  2  and 


ini  {T^-T^+L-l-^Y^/a-l )  ] 

ln[ dj-T^-fL-l )  (T<t-T>-«,-1-Q)/{t-1 )  (i2-r>+L-1-fi)] 

M, 

tn[  (2T2-T-VfB-«.-1  >/(L-1 )] 
tn[  (T^-T^+L-l  )/<L-1 )] 

tn[ (T.-T  +L) (T  -T  +L-1+ST  )/L(T.-T  +L-1)] 

I  •  •  I  O*  1  | 

ln[  (T1-»>-«.)  (Ij-T^+t-l )  (T--T>+L-1-Q)A(T1  ^T^+L-1  >  (Tj-T^+L-I-Q)  ] 

"« 

tn[  (2T2-'r„-*B*L-1 )  )] 

tnt  (T— )  ] 

tn[ 

"ll 

in[(2T1-'r--T>+L)/L] 
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now  occur  at  th« 


0(1)  location,  and  the  problem  btcoau  analytically  intractable,  in 


general,  if  they  arc  indistinguishable.  TO  circunvent  this  difficulty,  a  nodal  of  tha 
amrged  flames  is  oonsidarad  in  which  tha  two  reaction  sonas  ara  distinct,  with  thicknesses 
of  disparate  orders.  Tha  thinner  f laaa  is  embedded  in  tha  broader  flan*  and  appears  as  a 
discontinuity  on  the  seals  of  the  thicker  sons,  as  shown  in  Figure  3b.  This  nodal  was  also 
used  by  Xapila  and  Lodford  (11].  In  the  present  work,  the  f lanes  for  which  the 
deconposit ion-reaction  sons  is  the  thinner  one  are  referred  to  as  Type  A  f lanes;  those  with 
a  comparatively  narrow  diffusion  flans  are  called  Type  B  f lanes. 

A  third  kind  of  merged  flans  is  also  oonsidarad.  In  this  Type  C  flans,  both  reactions 
are  spread  over  the  sans  o(  1 )  interval  at  r, ,  but  the  structure  problem  becomes 
tractable  by  assuming  that  tha  right  and  left  side  of  the  fuel  equation  (2.1b)  are  of 
disparate  orders.  This  requires  that  Ij  be  snail  and  that  the  right  side  dominate 
aquation  (2.1b).  The  situation  thus  described  is  a  decomposition  reaction  in  which  very 
little  fuel  is  produced  relative  to  the  amount  of  nonopropellant  present,  and  whatever  fuel 
is  produced  is  immediately  consumed  by  the  bipropellant  reaction. 

The  kinds  of  structure  problem  which  result  from  the  models  described  above  are 
illustrated  in  Appendices  2  and  3,  which  consider  the  structures  for  solution  category 
15.  The  outer  solution  and  structure  results  for  that  category  are  as  follows. 

Assume  that  the  monopropellant  and  bipropellant  reactions  both  occur  at  r  ■  r,, 
with  1  <  r,  <  -,  and  that  T^r.)  -  o(1>,  Tf(r.)  -  o(1),  TQ(r,)  -  0(1).  The  0(1) 
solution  of  (3.2)  which  satisfies  these  assumptions  and  conditions  (2.2)  -  (2.5)  is  then 

1  <  r  <  r*t  Tm  -  1-e"/r*e"M/r,  T  -  rm~h*U,\~m/v  ,  (3.6) 

H|/r 

r,  <  r  <  -i  Jm  -  0,  T  -  T  -*»1«Qt(T--«sVfc-1^>e  .  (3.7) 


& 

s 


In  order  for  the  decomposition  reaction  to  take  place  at  r„  the  temperature  there 
must  be  Tj,  with  Tj  >  tf.  It  follows  that 


H/r.  V*  ♦L-1-Q 


T  -T  ♦L-l-Q 


H  fTlVl"/r* 

,  a  -  ( - r*— )e 


(3.8) 


1 

*1 


Consistency  with  the  hypothesis  end  nonnegative  ness  fractions  requires  either 

t  >  0  end  aex(T  ,T  -L+1+0.T  _)  <  T.  <  T 

9  I  B8  1  * 

or 

L  >  0  end  T,.  <  T,  <  *in(TB-L«-14Q,TBg)  . 

Since  Yg  cen  vanish  for  r  ^  r,,  the  bipropellent  reaction  can  occur  at  r,  if  t2 
does  not  exceed  T^. 

the  fleas  structures  aust  now  be  analysed  to  uncover  additional  constraints  on  the 


validity  of  the  above  results.  Appendices  2  a  3  treat  the  case  where  T2  is  less  than 
T1 ,  hence  only  Type  A  and  Type  C  flaaes  are  considered.  (Type  8  f lanes  do  not  exist  for 
T2  <  Tj,  as  aay  easily  be  shown.)  It  is  found  that  the  results  (3.6)  -  (3.8)  are  valid 


t  >  0»  T2  <  T1>  and  aax(Ta,TB-I.+  1+Q,TM)  <  <  Tm  or 


ttT.fT^-I.+O-K))/*)  <  T,  <  ain(Tg-L+1+fi,TBS) 


An  approach  similar  to  that  just  described  yields  the  results  shown  in  Table  3  for 
aerged-flaae  categories  14-17.  Details  of  the  analysis  are  given  by  McConnaughey  [10]. 


fable  3, 


M14  -  inKT^-T^-l-O/a-l-Q)] 

M15  -  InKT^^XT^-T^+L-l-iD/lCT^^+L-l-Q)] 
"16  ”  M13 

M17  -  InKT^-T^-tUAJ 


3.4.  H  »  1 

For  M  >>  1,  relations  (2.4)  and  (2.5)  become 

♦  Yf  -  Yc  -  1,  Ym  ♦  QY0  +  T  -  T,-L+1  ,  (3.10) 

for  r  <  “  while  equations  (3.2)  imply  that  YB  and  T  are  constant  outside  of  the 
f lanes,  which  is  incompatible  with  the  boundary  conditions  (2.2).  There  must  therefore  be 
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a  reaction  none  at  r  -  1.  In  addition,  aquation  (3.10a)  can  only  be  true  if  YQ(r>  “  0 

everywhere,  ainoe  the  sum  +  Yf  ♦  Yq  cannot  exceed  unity  and  Yq  >  0  (by  definition  of 

aaaa  fraction).  But  thin  violates  the  condition  Y  »  Y  at  r  -  «,  so  that  there  aunt 

o  om 

also  be  a  layer  at  r  -  •.  Clearly  the  reaote  reaction  auat  be  the  bipropellant  reaction, 
while  the  deccaposition  reaction  occurs  in  the  surface  layer  and  entirely  consuaes  the 
aooopropellant  there.  Hence, 

rm(ri  -  0,  T(r)  -  Tg-L+1  for  1  <  r  <  •  . 


The  analysis  of  tha  reaction  sonea  is  given  in  Appendix  4  and  leads  to 
M  -  *18  -  a^1(Tt-L-t1)2exp[ei(TB-l+1-T1)/2T1(T<-L-*-1)l 

which  is  valid  for 


0  <  L  <  1  and  T,  <  T^-L+1  . 


(3.11) 

(3.12) 


Considering  the  liait  of  large  activation  energies  yislds  an  0(1)  burning  rate  M 
in  solution  categories  1-17.  The  diaensional  evaporation  rate  N  thus  varies  with  the 
diaensional  droplet  radius  a,  i.e. 

M  -  ha  ,  (4.1) 

where  k  is  a  constant  with  appropriate  units.  Solution  category  18,  on  the  other  hand, 
gives  N  ~  for  sufficiently  large  1^'  in  which  case 

M  -  ka  2  .  (4.2) 

Mote  that  this  is  only  possible  for  L  <  1,  i.e.  when  the  heat  of  deccaposition  exceeds 
that  of  vaporisation,  but  in  practice  this  is  alnost  invariably  the  case  of  interest 
(Ludford,  Yannitell  and  Buckaaster  [12]). 

It  is  well  established  that  (4.1)  holds  for  pure  bipropeliant  burning  and  for 
deccaposition  of  saall  droplets  in  an  inert  ataosphere,  while  (4.2)  holds  for  simple 
decompositions!  burning  of  relatively  large  drops  (Williams  [13]).  experimental  evidence 
of  the  relation  between  M  and  a  for  the  hybrid  combustion  problem  considered  in  this 


work,  however,  is  limited  and  inconclusive.  Investigations  of  the  decoapoaitionsl  burning 
of  a  hydrasine  droplet  la  oxygen  do  auggeat,  though,  that  (4.2)  ia  aatiafied  for  all  but 
the  tlloet  droplets,  when  (4.1)  holda  (Dykema  and  Greene  (3]>  Roaaer  aad  Peakin  (5) » 
Allison  and  Faeth  [6]). 

An  iaoreaae  in  droplet  diameter  reaulta  in  an  inereaae  in  the  reaction  ratea  D1 
aad  Dj.  (Recall  that  it  ia  the  variation  of  M  with  and  02  that  ia  sought 
here.)  Definition  (3.1)  indicates  that  the  response  M(D^ ,  Dj )  may  be  described 
qualitatively  by  the  behavior  of  MlT^,?^).  The  latter  is  deduced  from  the  results 
reported  in  (3.11),  Table  2  and  Table  3  by  allowing  Tj  and  tj  to  vary  in  the  burning- 
rate  formulas  while  all  other  parameters  are  held  fixed. 

Different  responses  Wf^Tj1)  are  found  in  each  region  of  the  parameter  plane  shown 
in  Figure  4.  The  domain  of  the  response  associated  with  selected  regions  is  shown  in 
Figures  5.  Figures  5a-Sc  show  that  IKt'1,*”1 )  is  a  single-valued  function  when  Tw-T 
ia  positive  and  greater  than  -L-M+Q.  The  resulting  response  surfaces  are  found  to  be 

w  1 

continuous  and  monotonically  increasing  with  T^  and  ,  as  illustrated  in  Figures  6. 
Far  any  point  on  these  surfaces,  the  flams  temperatures  do  not  exceed  the  ambient 
temperature  Tm  and  heat  is  gained  from  the  atmosphere!  the  decomposition  at  r,  is 
couplets  and  all  fuel  is  consumed  at  r2. 

When  0  <  Tm  -  <  -t  ♦  1  ♦  Q  is  satisfied,  the  function  MtT^.Tj1)  is  multiple- 

valued  over  part  of  the  domain,  as  illustrated  by  Figures  5d-5f.  The  flame  temperatures 
are  no  longer  bounded  above  by  T„  and  heat  ia  now  lost  to  the  ambient.  Decomposition 
at  rj  is  complete  in  Regions  IV-VII »  it  is  only  partial  in  Regions  VZIZ-XIV  for  certain 
values  of  Tj  and  T2.  In  Regions  ZV  and  V,  the  fuel  is  entirely  oxidised  at  r2,  while 
portions  of  the  responses  associated  with  Regions  VZ-XZV  correspond  to  partial  fuel-burning 
diffusion  flames. 

The  nultivaleedneas  of  H(T,) 1 .Tj1 )  is  manifested  by  the  appearance  of  folds  in  the 
response  surface,  as  seen  in  Figures  7,  and  corresponds  to  the  phenomena  of  auto-ignition 
and  auto-extinction  of  one  or  both  reactions.  The  lower  turnaround  point  in  a  constant-^ 
(or  T2)  cross-section  of  the  surface  is  the  ignition  point  at  which  the  solution  jumps 
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Figure  Sc 


(See  caption  on  page  19.) 


from  a  nearly  axtinguishsd  it»t«  up  to  in  ignited  etete  ee  T~  (or  T1  )  ie  increased 
through  its  local  maximum.  The  upper  turnaround  point  ie  the  extinction  point  where  the 
solution  juspe  from  an  ignited  state  to  a  weakly  reactive,  nearly  extinguished  etete  as 
Tj1  (or  t”1 >  is  decreased  through  its  local  sininun.  The  siddle  level,  which  is 
bypassed,  is  believed  to  be  unstable  since  it  inplies  a  decreaaing  burning  rate  with 
increasing  reaction  rate  or  pressure  (Buckmastar  and  Ludford  (1]i  Kapila  [14]). 

The  responses  in  Regions  IV,  VI,  VIII  and  X  differ  only  in  their  middle  level,  hence 
the  realisable  responses  in  these  regions  are  the  same,  similarly,  Regime  V,  VI I,  IX  and 
XI  share  responsea  which  are  effectively  the  same,  as  do  Regions  xii-xiv.  Figures  7a  end 
7b  thus  indicate  for  Regions  IV-XI  the  possibility  of  ignition  and  extinction  of  both 
reactions  simultaneously  or  of  the  bipropellant  reaction  alone.  Figure  7c  shows  that  where 
ignition  and  extinction  occur  in  the  response  of  Regions  XII-XIV,  they  occur  in  both 
reactions,  ignition  being  possible  only  when  the  diffusion  flams  is  remote. 

4.2.  Comparison  with  Related  Works 

Three  "edges"  of  the  response  surfaces  found  here  represent  previously  addressed 
problems,  so  that  related  results  nay  be  compared. 

The  bipropellant  reaction  is  found  to  be  remote  for  0  <  t”1  <  min(T~\T~1)>  the 

2  BS  ® 

burning  rate  is  then  controlled  by  the  decomposition  reaction  alone.  The  corresponding 

edge  of  the  response  surface  therefore  describes  pure  nonopropellant  decomposition.  That 

problem  has  been  considered  by  Li  nan  [15]  and  by  Ludford,  Tannitell  and  Bucknaster  [12]  in 

the  limit  of  large  activation  energy.  For  Ta  <  Tm,  the  work  of  Ludford  et  al.  produces 

the  burning  rates  Mj,  H10,  H^,  M17  and  in  the  a  asm  regions  of  the  Tm  -  vs. 

L  parameter  plane  and  with  the  same  domains  as  found  in  the  present  work  for 

Tj1  <  min(TBg,T<>1 ).  similarly,  there  is  agreamsnt  between  the  associated  response  curves 

N(0^>  obtained  by  Ludford  et  al.  [12]  and  the  constant-Tj  cross-sections  of  the  response 

-1  -1 

surfaces  found  in  this  work  for  T.  <  T  .  Linen  [15]  obtains  the  above  expressions  for 


N,  but  does  not  specify  the  parameter  values  for  which  each  is  acceptable 


Hut,  it  My  be  noted  that  for  T,  <  and  L  >  1,  the  decomposition  flaae  is 
sitting  on  the  droplet  surface  and  the  problem  in  1  <  r  <  •  is  that  of  an  ordinary  fuel 
drop  at  temperature  T>f  but  with  latent  heat  L-1.  Appropriately  modifying  the  results 
of  law  [16],  obtained  Tie  activation-energy  asymptotics,  yields  burning  rates  N,,  Mj,  *4 
and  Nji  ths  corresponding  response  curves  N(D2>  are  consistent  with  the  constant-*^ 
cross-sections  of  the  response  surfaces  found  hare  for  L  >  1  and  <  Tg. 

Finally,  the  results  of  Buckmaster ,  Kapila  and  Ludford  [9]  are  considered.  As 
mentioned  earlier,  these  authors  used  DamkBhler-number  asymptotics  in  their  treatment. 
They  were  thus  abls  to  determins  only  ths  three  "corners" t 


lim  H(D  ,D  >,  lie  Him  M(D  ,D  >],  lim  llim  M(D  ,D  )] 
D1,D2*«  '  02*0  01**  D^O  D2+» 


of  the  response  surface.  As  expressions  ,  Nj,  and  k17  were  obtained.  Formula 

Mjg  was  also  found  by  considering  the  limit  ^  •  for  M  »  1.  More  precisely,  they 


M1S  for  l.  <  1 


lim  m(d1,D2)  -  /  M,  for  1  <  t  <  1  ♦  Q  ♦ 


mu  for 


r«18  **  *•  <  1 

Lim  tlim  M<D.,D2)]  -  < 

92*0  D,-  \  Mj  for  b  >  1 


lim  (lim  KID-.D.)]  -  M  _ 

V°  D2"« 


in  agreement  with  the  corresponding  corners  of  the  responses  found  in  the  present  work. 
Discrepancies  exist,  however,  between  the  nature  of  the  curves  connecting  these  corners  as 
conjectured  by  Buckmaster  at  al.  and  the  curves  dictated  by  the  finding  of  the  present 
work,  as  illustrated  in  Figures  8. 


lit*  problau  within  the  decomposition  flana,  given  by  (2.1n,d)f  (3.3)  and  (3.4)  in 
i2*m/4t2  -  -d2*t/dC2  -  xaD1«xp  xt  , 

dx^/dC  ♦  -*/r2,  dxfc/dC  ♦  M(T1*T-+l*)/r2  ns  C  ♦  —  , 

dx^/dC  ♦  0,  dxfc/dC  +  Mdj-T^+L-D/r2  at  C++*  , 

what* 

*?/«i  « ** 5i - e2,); - o(i>  ' 

»a  -  exa(C>  ♦  o(e),  i  -  ♦  ext(C>  ♦  o(e>  . 

This  aay  b*  rwwrittan  as 

2d2x_  dx  dxa 

- j-  -  xa*xp(«y-xa),  ♦  -1  as  y  ♦  -•#  ^r“  +  0  as  y  ♦  +•  ,  (Al.1) 

dy  y 

whar*  y  -  *C/r2  ♦  a"1[a0+tn(2D1r*/*2)],  a  -  (T^^+L-1)  and  ag  is  an  undetermined 
constant.  Boundary-value  problem  (Al.1)  has  baan  traatad  numerically  by  Linan  [17]  whoa* 
r* salts  indicat*  that  a  solution  exists  provided 

T1  >  Tg  -  t  ♦  1/2  if  x^  ♦  c  >  0  as  C++*  . 

*,>*,-  L  ♦  1  if  x^+0  as  C++*  • 

Diffusion-Flan*  gone 

Tha  following  bolds  in  the  diffusion  flan*. 

T  -  6b„  +  o(4),  T  -  T  +  6y  (n>  ♦  o(4)  , 

m  0  BS  t 

d2yt/dn2  -  -«OO^TfTo*xp(82A2-«2ABfi>  '  (*1>2) 

yt  -  bjtv+bj  as  h  ♦  — ,  yt  ~  b^-t-^  as  n  ♦  +•  , 

where  n  “  (r-r2>/6,  4  «  1,  bQ  is  a  nonnagativa  constant.  <  Tgg,  b1  • 
h(T--T-^>1-«)*xp(-«/r2>/r2,  b2  -  MtT^+L-1  )exp< -*/r2 >/r* ,  b3  and  b4  ar* 


undetermined  oonetente ,  end  Yf  end  Yg  ere  given  by  (2.4)  end  (2.5).  Coneider 
*2  «  T»S*  Then  ( A 1 . 2 )  any  be  written  ee 

d2y/dx2  ■  y  ~  i*  u  x  *  +•  ,  (*1.3) 


where 

y(n)  -  yt(n)  -  (bf+bj^-toj+bj)^,  x  -  (bj-b,  )n/2+(b3-b4>/2  , 

and  <diere  the  eeell  per  erne  ter  5  he*  been  ehoaen  according  to 

«3  -  l(b2-b1)/2)2(Dj)"1Q  •*P(92/»M-»2A2)  . 

Physical  ooneiderations  require 

y  <  -|x|  for  ell  x  .  (A1.4) 

equation  end  boundary  condition*  (*1.3)  constitute  the  standard  structure  problem  for 
a  Burke-Sehuaenn  flea*  which  has  been  shown  by  Bolaes  [18]  to  possess  exactly  two 
solutions,  only  on*  of  which  satisfies  (*1.4). 


Appendix  2i  Structure  Analysis  for  Category  15i  Two*  *  Flames 
Diffusion-Fleas  gone 

The  appropriate  expansions  in  the  bipropellant-reaction  son*  include 
r  -  r#  ♦  «n,  -  «ya(n>  ♦  o(«>,  T  -  T1*«yt(n)+o(«)  , 

where  i  2  ■  D'Y^exp^j/Tj-Oj/Tj),  To#  “  Yo(r,>,  »  expfcflj)  with  c  >  0,  yfc  <  0  for 

n  <  n#,  y t  "  0  at  n#,  yt  <  0  or  Y^  ■  0  for  n  >  n#  with  n*  the  location  of  the 
deooapoeition  reaction.  The  function  ya(n)  thus  satisfies 

d2y  dy  dy 

— ~  -  o(  i )  for  n  J  n.,  “  ♦  -  **  n  ♦  ♦  o  as  n  ♦  +*  . 

>  w  dn  _2  dn 

on  r# 

Continuity  of  ym,  yt  and  d(ya+yt)/dT\  at  na  is  required,  therefor* 

2  2 
ya  *  c^-nn/r*  for  n  <  ya  -  Cg-un./r,  for  n  >  n#  , 

yt(%)  -  yt(n.>,  3*  (n  )  -  -*  (n.)  ♦  -j  , 

re 

where  the  constant  ag  is  undetermined.  The  equation  for  yt  bscoaes 


(A2.1) 


df yt/dn  ■  yt-e1n-c2  for  n  <  n,#  d  yfc/dn  -  yt”c3l'~c4  for  n  >  n*  *  (*2.2) 

whara  c1  -  iKT^^+D/r*,  c3  “  Ht^-T^+L-l-O/rf ,  -  (Cj-Cjln.+Cj  and  c2  ia 

dataminad  by  sntching  higher— ordar  taraa.  Iha  solution  of  (A2.2)  mat  satisfy  (A2.1)  as 

wall  as 

dyt/dn  ♦  c1  as  n  ♦  dyt/dn  ♦  e3  as  n  +  . 

It  follows  that 

n-%  n#-n 

yt  ■  CjS  ♦  c^n  ♦  c2  for  n  <  yfc  -  c5a  ♦  e3n  +  e4  for  n  >  n*  , 

2 

whara  cj  -  -QM/2r# .  Iha  raquiraaant  y  <  0  for  n  <  is  tharafora  aat  provided 
*1  *  Ts  "  1  +  Q/2‘ 


Deocnposttion-Flane  Iona 

In  tha  aooopropallant  raaction  at  H#,  tha  following  is  valid. 

r  -  r.+«n.*«e,  -  dy^  +...♦  cx^CO+oU),  T  ■  ♦  cx^C)  ♦  o(e)  , 

whara  c  «  yB(n.),  and  tha  eaaponants  of  YB  pracading  tx^  ara  constants 

For  y.^  f  0,  tha  problan  for  is 

d2x  x  dx  dx 

“  “V  *  C3**1  *  0  “  5  *  d^  *  "C5  *  C3  <  0  "  *  *  *”  ' 

#w 

whara  ■  cdy^  D*  ■  0(1).  Zntagration  shows  that  Cj  ”  *Oj. 

For  ”  o(0  at  (|(#  tha  structure  prablasi  is  aqoivalant  to  (*1.1)  with 
«  -  T1-Y#*L-1-Q/2. 

It  follows  that  a  Typa  *  nargad  flans  nay  axist  for  catagory  15  if  0^  »  expCcO^) 
with  c  >  0,  and  if  nax(Ts-L+1*Q/2,T2 )  <  T,,  <  T„  or  naxlT^.T  -L+(14Q)/2,T2 J  <  *r 


Ths  problem  in  ths  Typ#  C  flat  conn  for  category  15  is 

d^x  dx  dx 

■— j5  -  O1x>«cptc35*d0-(1'*Q>x|||]<  ♦  -  **j  *■  5  ♦  "*»  *  0  **  C  ♦  , 

dT  r* 

where  -  cx^tC)  ♦  o(e),  *  ■  *3  ♦  ext<5)  ♦  o(e),  5  -  (r-r„)/e,  e  -  T2^  * 

01  -  c  O'  *  0(1)  and  dg  i*  an  undeterained  constant.  This  nay  be  rewritten  in  the  fora 

of  (A1.1)  and  a  solution  thus  shown  to  axist  provided  T^  >  aax [T^.T^-L+f 1+Q)/2J  or 

T^-L+1-KJ  <  T ,  <  I„. 

The  saall  Mount  of  fuel  produced  in  the  decomposition  reaction  is 
Yf  -  eD'^x^e  VlOjI^ejqpCS^Aj-^jA^  >1  /  which  is  indeed  o(e)  since  T^  >  Tj  and 
»  1. 


Appendix  4t  Structure  Analysis  for  Cat 


Surface  Laver 


The  region  adjacent  to  the  droplet  is  described  by  C  ■  y  (r-1)  with  y  <<  1,  and 
the  burning  rate  M  is  of  the  fora  M  *  Y  '*  where  H  “  0(1).  Within  this  layer,  the 
following  holds. 

Y  -  a(U  ♦  o(1),  T  -  t(C>  +  0(1),  t  -  T-L+1-a  , 

a  S 

2 

-  W  £  ■  r2D'a  exp(8  /f  -8  /t)»  *  H(a-1)  -  -HI,  a  -  1-L  at  (*0  . 

«  1  ill 

Let  C.  label  the  saallest  value  of  C  at  which  t(C#)  ”  t#  -  aax(t(01»  take  y  “ 
(D'exp(81/r1-81/t#)]  and  require  T^  <  t*.  Than  the  deccaposition  reaction  occurs  at 

r  •  1  ♦  YC,  and 

a  -  1-exp (H(C-C#)]  for  ?  <  C#,  a  *•  0  for  C  >  C,,  C*  *  “*  ^nL  . 

The  structure  problem  at  is  given  by 

dJx 

— r®  ■  D  x  exp(-x  ),  x  ~  -HC  as  C  ♦  -•»  x :  ~  0  as  5  ♦  +•  , 

»  ■  a  a  a 

where  C  -  U-C.)/e,  e  *  t*/8  ,  T  -  ex .(C)  +  o(e),  T  »  t#-ex  (5)  +  o(e)  and 

w  »  a  a  a 


*-  JU  JV  w>  -NlV-Vn  «  e^  a  «  a  O*  •  «> _s  «  m.  CaVjs  "-a.  'a  a  ■  ’-«■  a»^A - 


D1  -  c  D|  -  0(1).  Thar  afar*,  N  -  /2D1  which  specifies  N.  This  result  is  valid 
provided  L  <  1  sad  T1  <  Tg  *  L  ♦  1. 

Tha  ooordinats  of  the  reaote  ccabustion  field  is  *  -  Tr.  sad  the  bipropsllaat 
reaction  is  located  at  X..  The  tsaperature  equation,  given  by 

d2T/d*2  ♦  ( ( 2im)/R2]  dT/dR  -  -Y"2QDjTfToexp(e2/»2-«a/*>  . 
is  reactionless  for  *  »*  «*  and  for  say  »  1,  Oj  »  1  provided  the  product  Tf»0 
vanishes  for  It  /  X..  The  bipropellant-reaction  zone  is  consequently  a  Ber* e-Schwann 
diffusion  fleas  with  TQ  -  0  for  *  <  *,  and  Yf  -  0  for  *  >  Continuity  at  S. 

fixes  the  value  of  X.. 

The  structure  problea  at  R*  has  the  fora  of  (At. 3),  hence  a  solution  exists  sad 
additional  restrictions  need  be  lapoeed. 


List  of  Symbols 


dimensionless  quantitiss : 


mss  fraction  of  species  A;  A  -  a,  f  or  o 
tesiperature 


_  2  i 

DamkShler  number,  varies  with  a  p  >  i  -  1  or  2 


activation  energy 


ratio  of  the  heat  of  combustion  of  the  bipropellant  reaction  to 
that  of  decomposition 


Ti 


evaporation  rate  (“burning  rate") 
radial  coordinate 
latent  heat  of  evaporation 
»  Da  expC-O^/r^ 

positive  parameter  characterising  the  magnitude  of 


dimensional  quantities i 
a  radius  of  droplet 

p  pressure 

M  evaporation  rate,  varies  with  all 


subscripts  s 


monopropellant 

fuel 

oxidant 

pertaining  to  decomposition  reaction 
pertaining  to  bipropellant  reaction 
surface  value 


ambient  value 
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